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Abstract 

The anomalies observed in recent cosmic ray experiments seem to strongly con- 
strain the nature of the dark matter. In this letter, we investigate a possibility of the 
fermionic dark matter with a minimal extension of the standard model. We found 
that the dark matter decays caused by the dimension six operators can naturally 
explain the anomalies. 



1 Introduction 



The presence of dark matter (DM) has been established by numerous observations, which 
requires physics beyond the standard model (SM). The nature of the DM is now one of 
the most important issues not only in cosmology but also in particle physics. 

The minimal extension of the SM is to introduce one extra particle called X, that is, a 
candidate particle for the DM pQ. We assume the DM X to be a singlet of the SM gauge 
groups. The X particle can be a boson or a fermion. We consider in this letter the case 
of fermion. We denote the DM X as $ to distinguish it from the case of boson 

The new particle \I/ may have Yukawa couplings as 

X^tH + h.c, (1) 

which induces too fast decay of the Here, £ is the lepton doublet, H the Higgs doublet 
and i = 1,2,3 denotes family indices of the leptons. Thus, we assume the Yukawa coupling 
is strongly suppressed, otherwise the \l/ can not be a candidate of the DM. The required 
suppression may be easily obtained by suitable configurations of particle wavefunctions 
in a higher dimensional theory 

If the above dangerous Yukawa coupling constants A; in Eq. (TjQ) are strongly suppressed 
as | Ai| <C 10~ 26 , the dominant operators become dimension six four- Fermi interactions 
among the SM particles and 

^We\ ^JW, (2) 

where i, j, k = 1, 2, 3 and M* is the cut-off scale, which would be regarded as the Grand 
Unification scale. The lifetime of the \l/ is given by 

r*~ 10 26 sec —r, * - , 3 

VlO 15 GeVV V / 

where is the mass of the DM We see that the lifetime is much longer than the 

age of the universe in a large parameter space of M* and m# and hence the \I/ can be a 

candidate of the DM. 



1 See Refs. [2] for the case of X being a hidden gauge boson A' . 

2 An example is given in a 5-dimensional space-time with an extra dimension S 1 /Z?. We put the ^ 
on one boundary and the Higgs H on the other boundary. We put quarks and leptons in the bulk. Then, 
we find the Yukawa coupling ^!lH is exponentially suppressed if the size of the extra dimension is large 
enough, compared with the inverse of the cut-off scale in the 5-dimensional theory. 
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The \l/ may be produced non-thermally in the early universe through the above di- 
mension six interactions. The density parameter of the ^ is given by 

0*' = 0(0.1) (-^M* ( 10 ' 5 GeV V (J5S-V (4) 
^ V ; UO n GeW V M * J TeV ^ 

where Tr is the reheating temperature after the inflation We see that the observed DM 
density f^DM^ 2 — 0.1 is also explained for a wide region of the parameter space of M*, m$ 
and Tr. 

The purpose of this letter is to show that the anomalous excess of cosmic ray electron 
and positron recently observed by PAMELA [3] and ATIC [Ij/PPB-BETS [5] is naturally 
explained by the decay of the DM \I/ Q In the present analysis we assume the following 
dimension six operator, for simplicity; 

j^m^f) + a (v^)(W)] + h.c, (5) 

* 

where a is a free parameter of 0(1). We find that the replacement of the third family by 
the second family does not significantly change our final result and hence we consider only 
the above operators. However, the more general analysis including all possible dimension 
six operators is straightforward and will be given elsewhere. 

For completeness, we also study the case that the two-body decays of the DM are 
induced by Eq. ([I]). We find that the required lifetime for the DM \1/ is given if the 
Yukawa couplings |Aj| are sufficiently small as (9(10~ 26 ). We choose the flavor structure 
as | Ai| |A 2i3 | in the present analysis. 

2 Cosmic rays from the hidden fermion DM decays 

Let us discuss the cosmic ray signals from the decays of the DM The interaction ([1]) 
dominantly causes two-body decays of the DM 

^ -> hu, Zu, W ± e T (6) 

with branching ratio 1:1:2. 



The derivation of Eq. ((4} is as follows: The production cross section of the ^ is (av) w T 2 M^ 4 . 



Following the Boltzmann equation, one can get /n ra d ~ n ra d (<rv) H 1 



t=t r 



'For recent progress in the study of the decaying DM signal, see Refs. [U |H1 H] • 
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The interaction (jSJ) mainly causes three-body decays 

* -> r ± e T z/, e ± e =F z/, dfru, J&u (7) 

with branching ratio 1:1: 3a 2 : 3a 2 . Here, we have assumed that the \I> is a Majorana 
fermion and that ~ 110 GeV and my 3> m^. In both cases Eqs. (EI) and (I7j), high 
energy electron, positron, photon and antiproton are produced. To estimate the energy 
spectrum of these decay products, we have used the program PYTHIA [10] . The particles 
produced in the DM decays are influenced by various factors in the propagation. For the 
propagation in the Galaxy, we adopt the method discussed in Refs. [H] with Navarro, 
Frenk and White halo profile [T2] ; 



Pdm — 7—7 — Ul . / , — ^77, [o 
[r/r c )[l + (r/r c )J 2 

where po = 0.26 GeVcm -3 and r c = 20 kpc. 
Positrons and electrons 

As a diffusion model, we use MED model in Ref. [13]. We see that the positrons and 
electrons come from the DM decays inside the Galaxy, and especially the decays within 
a few kpc from us are important. In Fig. [TJ we show the total flux of the electron and 
positron. The left figure shows the case that the interaction (jSJ) is a dominant interaction 
causing three-body decay. We set the DM mass = 1800 GeV, the lifetime = 9 x 10 25 
sec and a = (2a/2) _1 . Hereafter we fix a = {2^/Tj- 1 . The lifetime is given by ~ 3 x 10 
GeV in Eq. ([3]). For the right figure, the case that the interaction ([1]) is dominant is 
shown. Here, we set the DM mass = 1200 GeV, the lifetime r$ = 8 x 10 25 sec. The 
lifetime is obtained by Ai — 1 x 10~ 26 in Eq. (JT|). As for the background flux, we set 
0.0253(^/1 GeV)" 3 - 206 GeV^cm-V^r" 1 in both cases 3 . 

Next, we estimate the positron fraction. For the background, we extrapolate the 
background estimated in Fig. [1] and assume that the background consists entirely of the 
electron, since the secondary positron would be negligible for E ^ 10 GeV. Fig. [2] shows 



5 This background is estimated by fitting the data points as BG + signal, assuming BG is power-low 
and adopting the signal in the case of three-body decay. As the fitting parameters, we have used the 
lifetime of the a in Eq. ([5]), and the coefficient and power of the background. Here, we set the weight 
for each data point 1. 
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Figure 1: Positron and electron fluxes with experimental data [U [5]. Left: three-body 
decay. Right: two-body decay. The solid line represents the DM signal plus background 
and the dashed line the background. 
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Figure 2: Positron fractions with experimental data [2l[T3l[T5]. Left: three-body decay. 
Right: two-body decay. 



the positron fraction. One can see good agreements with the experimental data in both 
cases except in a low-energy region. The behavior of the background in the low-energy 
region is complicated due to various factors such as solar modulation or contamination 
from the secondary positron. This can be reasons why the naive background estimation 
is not good in the low energy region. The detailed treatment of the background is out of 
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the reach of this letter. 



Gamma ray 

For the gamma ray, both of the DM decays in the halo and extra-Galaxy are important. 
To estimate the halo component, we have used the NFW profile in Eq. (JHJ) and averaged 
the halo signal over the whole sky excluding the region within ±10° around the Galactic 
plane. 

For the extragalactic component, the gamma ray is influenced by the red-shift. We 
estimate the extragalactic component by using the following cosmological parameters; 
tt*h 2 ~ 0.11, ft m atter/* 2 ^ 0.13, fi A ~ 0.74, p c ~ 1 .0537 x 10- 5 h 2 GeVcm" 3 , h ~ 0.72 [16]. 

In Fig. [3l the gamma ray fluxes are shown. We set the background flux as 5.18 x 
10~ 7 (E/1 GeV) -2 ' 499 GeV _1 cm~ 2 s _1 sr _1 as in Ref. [17]. We have assumed the energy 
resolution is 15%. In both cases, the DM signals are consistent with the current experiment 
data and anomalous behavior of the gamma ray flux is expected to continue up to higher 
energy. This will be tested by upcoming Fermi Gamma-Ray Space Telescope (FGST, 
formerly named GLAST |18j). 
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Figure 3: Gamma ray fluxes with experimental data [HJ [20]. Left: three-body decay. 
Right: two-body decay. Solid line represents the DM signal plus background and dashed 
the background only. 
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Antiproton 

We estimate the antiproton flux, following the Refs. [9]. For the solar modulation, we 
set 4>f — 500 MV. In Fig. HI we show the antiproton fluxes for some different diffusion 
models [13]. Here, we show only the DM signals. We can see the contradiction between 
the experiments and the signals in some diffusion models. However, in both cases of the 
two- and three-body decay, MIN models (and also MED model for the three-body decay 
case) do not conflict with the experimental data. Therefore, for the antiproton, the DM 
ip is consistent with the experimental data at least in some diffusion models. 
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Figure 4: Antiproton fluxes with experimental data [ZD 122] . Left: three-body decay. 
Right: two-body decay. Solid, dash-dotted and dashed line represent MIN, MED and 
MAX diffusion models, respectively. 



3 Discussion and conclusions 

In this letter, we have investigated the case of the fermionic DM. Both two- and three- 
body decays can explain the anomaly of electron and positron cosmic ray. In addition, 
they are also consistent with the gamma ray anomaly. Especially, the decay caused by 



the dimension six operators seems to be attractive, since it natural 
lifetime of the DM for the GUT-scale cut-off M, ~ 10 15 -10 16 GeV 6 



y explains the proper 



'DM decays via GUT-scale physics have been also discussed in a recent work 
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There are remaining issues. For example, the reason of large suppression of Eq. (CQ) is 
unclear. In addition, the reason why the DM ^ decays dominantly into the first family 
leptons is unclear. However this problem may be solved by choosing suitable Froggatt- 
Nielsen charge for the DM ^ 3- 

Finally, we should note that the hidden fermion \l/ can be identified as the lightest neu- 
tralino in the supersymmetric standard model. In this case, R-parity breaking operators 
may induce the DM $ decays in similar ways discussed in the present letter. 

Note Added: Discrimination between two- and three-body decay of the DM would be 
possible in the future cosmic ray experiments [26J. 
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